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Analysis of the singularities of the polarization field of CMB, where polarization is equal to zero, 

00 

0^ ' is presented. It is found that the classification of the singular points differs from the usual three 

types known in the ordinary differential equations. The new statistical properties of polarization 

C . 

^ ^ field are discussed, and new methods to detect the presence of primordial tensor perturbations are 



indicated. 
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I. INTRODUCTION 
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O ' electrons in cosmic plasma. The corresponding measurements of the CMB polarization are planning to be performed 



In the coming years the measurements of the angular anisotropics of the intensity of cosmic microwave background 
(CMB) by the cosmic missions MAP and PLANCK will possibly present one of the most promising methods of 
studying early universe as well as of precise measuring of basic cosmological parameters(see e.g. and references 
therein). In addition to the anisotropics of the intensity, it is possible, though more difficult, to measure polarization 
of the radiation. The polarization is a secondary effect induced by the scattering of anisotropic radiation field on 



in the coming space missions 



Polarization field is described by a traceless 2 x 2-matrix which can be decomposed into a sum of three Pauli 
. . matrices (Ji with the coefficients known as Stokes parameters: 



■ a = eifTi (1.1) 

^ : As is well known circular polarization does not arise in Thomson scattering (because of parity conservation), so that 
■ ■ ' ^2 = and the matrix a is symmetric. Usually it is parameterized in the form: 

The sources of polarization are anisotropics of radiation field induced by different types of perturbations, namely 
scalar, tensor and vector ones. Vector perturbations decay in the early universe but may arise at small scales at later 
stages and influence the CMB polarization in the case of reionization. 

Geometrical properties of the polarization field allow to obtain an important cosmological information [2-9]. An 
importance of the study of singular points of polarization field, where Q = U = 0, was emphasized in ref. pT| ]. In 
this paper geometric classification and statistics of the singular point was proposed. This has been done in terms of 
the fields Q and U which directly enter polarization matrix. The latter are 2-dimensional tensor fields and have the 
appropriate transformation properties under rotation of coordinate system. 
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II. TYPES OF SINGULAR POINTS. 



In this paper we will investigate the classification of singular points of eigenvectors of the polarization matrix a. 
Though the positions and statistics of the singular points remains the same, their types could be quite different. The 
eigenvalues are easily found: 



Ai,2 = ±VQ^+t72 (2.1) 

and the eigenvector corresponding to the positive A is 

ri(i) = (n,, ny) ^ (u, VQ' + ~ q) (2.2) 



This vector determines direction of maximum polarization and up to a normalization factor coincides with the direction 
of the vector P considered in refs. or orthogonal to it, depending upon the sign of the coefficient functions. 

The behavior of the vector field rV^^ in the vicinity of the singular points of polarization, Q = U = 0, is determined 
by the equation 



dx fix U 

Analysis of singular points of differential equations when both numerator and denominator can be expanded into 
Taylor series is well known and can be found e.g. in ref. fl^ j. In the usual case the following singular points can exist: 
focus, saddle, and knot. In our case the situation is more complicated due to the presence of the square root in the 
numerator which is generically non-analytic in the points where Q = U = Q. At this stage a question may arise why it 
is assumed that Q and U are analytic functions expandable into Taylor series (at least up to first order terms) around 
the points where Q = U = i), while the component of the eigenvectors are not. The reason for that is the following. 
The matrix elements of the polarization matrix Q and U are directly related to the anisotropy of the cosmic microwave 
radiation through the amplitude of photon-electron scattering. We do not expect any grounds for these quantities to 
have a square root singularity where their first derivative tends to infinity. On the other hand the eigenvectors of the 
matrix a just mathematically contain the square root Q"^ + U"^ and so it is singular dX Q = U = 0. The analysis 
of singularities of the vector field n can be done as follows. We assume that Q and U are expanded near singular 
points as 

Q = aix + 

U = bix + b^y (2.4) 



In the case when the matrix 

A/r — I 

^ bi 62 

is not degenerate, det M 7^ 0, it is convenient to introduce the new coordinates: 



M = K! 1: ) (2.5) 



x = Ai^ + A2r), y = Bi^ + B2rj (2.6) 

Evidently the types of the singular points do not change under this coordinate transformation. It is simpler to make 
the further analysis in polar coordinates: 

^ = pcos(j), 77 = psin0 (2.7) 
Equation ( ^.3[ ) in this new polar coordinates can be rewritten as 

where t = tan((/)/2) and 

N = -A2t^ + t\B2 - 2Ai) + t{2Bi + A2) - B2 (2.9) 
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and 



D = Ait^ - t^{Bi + 2A2) + t{2B2 - Ai) + Bi (2.10) 

Barring the degenerate case of Ai = we may take Ai = 1 without loss of generality. 

The behavior of singular points depends up on t he roots of denominator D. Let us first consider the case when it 
has three real root, ti,2,3- The solution of eq. ( |2.8D in this case can be written as: 



- - + 1) n u - <j p''^ (2.11) 
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where rg is an arbitrary constant and the powers Vj are 

= -i^±3%±l3h) (2.12) 

and so on by cyclic permutation of indices. Since it can be shown that = points where —> 00 are not 

generally singular. It can be easily checked that either all fj < or any two of them are negative and one is positive. 
In the first case the singular point resembles the usual saddle with the only difference that there are three and not 
four, as in the usual case, linear asymptotes/separatrices (see Fig. la). We will call it also "saddle". If one of i^j is 
positive (say vi > 0), and thus r becomes zero at t = ti, the behavior of the direction field near this point is quite 
different from the usual ones. The field line cannot be continued along (/) ~ (pi into = 0i + tt as can be done in the 
usual case. We will call this type of singularity a "beak" (see Fig. lb). 



In the case of one real root of denominator D the solution has the same form as (2.11) but now e.g. the powers 1^2 
and z/3 are complex conjugate. The solution can be written as 



- = (t^ + l) \ t-t2 1*^""^ exp(4/3Imi.2)(i-ii)'"' (2.13) 
ro 



where f3 = tan^^ [Imt2/(i — Ret2)]- The real root vi is negative, as is seen from eq. (2.12) and thus r does not vanish 
in vicinity of such singular point. The polarization direction field for this case is presented in Fig. Ic. This type of 
singularity can be called a " comet" . 

One can estimate the density of different singular points in the following way (see e.g. [| 11|). All singular points 
correspond to the case when both Q = and U = Q. The number density of these points is proportional to 

Qx Qy 



dQdU = det |^ j dxdy (2.14) 

and thus the density is given by the average value of the determinate, d = QxUy — QyUx- It can be easily checked that 
the saddle-type singularity takes place if d > that is the same condition as for normal saddles in the field determined 
by the equation dy/dx — Q/U see Q. It can be shown that saddles make 50% of all singular points (us) = 0.5{n), 
where n is the number density of all singular points. Calculations of the number of beaks and comets are more 
complicated and can be found numerically. According our estimates the surface densities for beaks and comets are 
correspondingly (ni,) « 0.052(n) and {uc) ~ 0.448(n}. We note that the probability of appearance of saddles, beaks 
and comets for random choice of Q^, Qy, Ux,Uy, is correspondingly Ws = 0.500, Wb ~ 0.116, Wc ~ 0.384. 

III. TYPES OF PERTURBATIONS FIELD AND TYPES OF SINGULAR POINTS. 

As we have already mentioned the polarization of cosmic microwave radiation arises due to anisotropy of the 
radiation field. It is a linear functional of the field and in the case of scalar perturbations the only way to construct 
two dimensional tensor quantity is to use second de rivat ives of a scalar function ^E*, as is argued e.g. in ref. §, [p^ . 



The matrix elements of traceless symmetric matrix (1.2) are constructed uniquely as 

a,y = 29^9^-^' - (3.1) 

For the functions Q and U it gives: 

U = 2dxdy^ (3.2) 



3 



In principle one may use also the invariant two-dimensional antisymmetric (pseudo)-tensor but parity conservation 
prevents from its appearance in polarization matrix in the case of scalar perturbations. For tensor perturbations there 
could b e sp ecific "external" directions in the spa ce a nd parity considerations do not prevent from using in the 



matrix (1.2) (see below). For this specific form (3.2) of polarization matrix there exist a particular (pscudo)scalar 
quantity which vanishes in the absence of gravitational perturbations |^,^ : 

B = eijdkdiMkj (3.3) 

Written explicitly it reads 

B={dl^ dl)U - 2d^dyQ (3.4) 

It evidently vanishes for Q and U given by expressions (^^). 

In the case when tensor perturbations are present, there is more freedom in polarization matrix a and the terms 
proportional to are permitted and the symmetric matrix a may be expressed through second derivatives of two 
independent scalar functions: 

My = 2d,dj-^ - S,,dH + e.kdkdj^ e,kdkd,<i> (3.5) 

This is a general decomposion of symmetric traceless tensor in 2 dimensions (see e.g. ^). 
With inclusion of tensor perturbations the functions Q and U take the form: 

U = 2d.,dy-^ - {dl ~ 9^)$ (3 6) 

The scalar B is expressed through the fourth derivatives of $ as following: 

B = 9^$ (3.7) 

The difference between scalar and tensor perturbations appears only in the fourth derivatives of the generating 
scalar functions $ and 5*, while the structure of their singularities is determined by the third derivatives. Thus the 
types of the singularities are the same for both types of perturbations. There are statements in the literature that in 
the case of scalar perturbations vector n cannot have curl, on the other hand tensor perturbations do produce a curl, 
see [|l3|, section 4. 

However the polarization tensor is not a vector but a second rank tensor and direct analogy with a vector field 
is not applicable. In the general case of singular points considered above the curl is not equal to zero for any type 
of perturbations. An explicit example of scalar generating function being a function only of r = i/a^+y^ near 
the singularity point shows that the latter may be either center or knot. These points are absent in our list of three 
presented above due to a specific degeneracy of this example. 

Thus to summarize, the singular points of the vector field n^^\ which is the eigenvector of polarization matrix 
corresponding to the direction of maximum polarization can be generically of the following three types (see above): 
saddle and beak (with three separatrices), and comet (one separatrix). In degenerate case, when some of the coefficients 
or their combinations (like e.g. determinants) may be zero, then there could be some other types singularities which 
we have not considered here. 



IV. STATISTICS OF SINGULAR POINTS. 

The functions Q and U are usually assumed to be independent Gaussian variables with equal dispersion Q : 

{Q')^{U')^al (4.1) 
Their first derivatives are also independent and non-correlated with the functions with dispersion 

{Q^QJ) = {U,UJ)=^,Jal/^ (4-2) 

where Qi — dQ/dx^, etc. All other correlators are zero. However the fact that in the case of scalar perturbations 
both functions Q and U as well as their derivatives are determined by a single generating scalar function ^' imposes 
some conditions on the correlators of the second derivatives. In particular the dispersions of the second derivatives 
Qij and Uij are not equal and these fields are correlated. The list of nontrivial correlators is the following: 
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{Qlx) = {Qly) = Y^0-2i {Uxx) = (Uyy) = Y^<^2i 
1 2 3 

{UxxUyy) = (UxyUxy) = — - 

{QyyUxy) ~ {QxyUyy) = ^^g- 
{QxxQ) = {QyyQ) = {UxxU) = {UyyU) = -all2, 



{QxxQyy) — {QxyQxy} {UxxUyy) {UxyUxy) i ^^2 ' (^■3) 

{QxxUxy} {QxyUxx} Tp^2' {QyyUxy} {QxyUyy} 1^*^27 



where erf = {{Qxx + QyyY) = {{Uxx + Uyy)"^) . Note that the asymmetric correlators in the third Hne are non- vanishing. 

These properties permit in principle to discriminate and detect a contribution of tensor (or vector) perturbations 
into polarization of CMB by measuring the dispersion of second derivatives of the Stokes parameters. In particular 
for pure scalar perturbations one should expect 



{Qxx) _ {Qly) 

{UD ~ {U^y) 



^ = 7/5, {{Uxx-Uyyf)=A{Ql\ (4.4) 



A deviation from this number would indicate a contribution from other, different from scalar, types of perturbations. 
The last equality in (27) corresponds to i? = 0, see Eq.(20). The property S = as a test for the absence of scalar 
perturbations was indicated in previous works, see for example Q. 

An interesting quantity which allows one to relate global characteristics of random field to local properties is the 
Euler characteristic, xe, see As it was noted in ref. this value is closely linked to the critical value of the 
amplitude of polarization, P = \/ Q^ + U^ , for which the regions of high polarization percolate. According to this 
[14] the percolation begins at the amplitude of polarization corresponding to xe = 0. 

This critical amplitude was estimated in ref. pT[ | where it was found that percolation takes place for p — 1, where 
p = P I oq is the dimensionless amplitude of polarization with unit dispersion. We estimated this quantity in somewhat 
different way than it was done in ref. [ pT| . In the 2D case the required value is defined by an equation 

Xe oc {pxx + Pyy)pe^v{-p^ /2) oc [p^ - 1) cxp(-p^/2) (4.5) 

that is identical to result obtained by ref. |^l|. Though the statistical distribution we used is different from the one 
of ref. [0, we got the same result: percolation occurs at p — 1. Let us remind that for the 2D Gaussian fields 
Xe oc pexp{—p^/2) and percolation occurs at p = 0. As it is following from (17) - (22) the same results are valid for 
all three types of perturbations of polarization field. This topics will be discussed in more detail in jl^ . 
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FIG. 1. Integral curves for three different types of singular points: 
peculiar solutions (separatrises) . 



saddle, b. 
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beak, & c. 



comet. Long dashed lines show 
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